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manifested when observables are related to mutually unbiased bases. In this paper, 
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essential feature of quantum mechanics, quantum entanglement. It is shown that a 
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of entanglement, independently of the choice of the set. This has implications for 
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basic properties of mutually unbiased bases composed only of product states. The 
constructions are illustrated with explicit examples in low dimensions. We believe that 
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also expect that they will be relevant to applications of bases in the experimental 
realization of quantum protocols in higher-dimensional Hilbert spaces. 
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1. Introduction 

Quantum complementarity forbids the simultaneous knowledge of almost all pairs of 
observables. This impossibility is drawn to the extreme in the case of observables 
described by operators whose eigenstates form mutually unbiased bases (MUBs). Two 
bases are said to be unbiased if any vector from one basis has an overlap with all vectors 
from the other basis that is equal in modulo. The definition for a bigger set of MUBs 
means that the unbiasedness property holds for all pairs of these bases. Accordingly, if 
we can perfectly predict a measurement result of one such observable corresponding to 
an eigenstate in one of the bases, then the results of all other observables corresponding 
to all other basis vectors of all other bases in the set remain completely uncertain. One 
typical example of a set of three MUBs is the eigenbases of spin-| projections onto three 
orthogonal directions: a spin-i state along one axis leaves us totally uncertain about 
the results along the orthogonal axes. 

A spin- 1 particle is a two- level quantum system, a qubit, and clearly admits three 
MUBs. A (i-level quantum system, a qudit with pure states described in d dimensional 
Hilbert space, can have at most d + 1 MUBs [1], and such a set is referred to as the 
complete set of MUBs. The first explicit construction of the complete sets of MUBs 
was presented by Ivanovic for d being a prime number [2]. Subsequently, Wootters and 
Fields constructed the complete sets for prime-power d [I]. Since then, many explicit 
constructions have been derived and they are collected in a recent review [3J . If is not 
a prime power, the number of MUBs remains unknown although it is considered unlikely 
that a complete set of MUBs exists in these cases. For example, the works [H El El [7] 
describe failed numerical attempts to find a complete set of MUBs in dimension 6. In 
addition to this fundamental question, MUBs find applications in quantum tomography 
[Ij, quantum cryptography [HI El [10], the Mean King problem [HI [121 [131 [H] , other 
tasks. 

Here we study the properties of entanglement between subsystems of a global system 
with a composite (i.e. nonprime) dimension as well as entanglement distribution among 
the states of MUBs. We show that the amount of entanglement, as measured as a 
function of the linear entropy of a subsystem, present in states of a complete set of 
MUBs of a composite dimension always must have a nonzero value that is independent 
of a chosen set. In other words, entanglement is always present in such a complete 
set of MUBs and it is always the same independent of the choice of the complete set, 
being solely a function of dimensions of subsystems. Moreover, for global dimensionality 
that is big enough, practically all MUBs of a complete set contain entanglement. We 
then show an experiment-friendly procedure that creates complete sets of MUBs in all 
dimensions d = p^, which are squares of a prime number. This procedure uses only one 
entangling operation, which is repeatedly applied to states of product MUBs to give 
the complete set. Remarkably, the generated set consists of either product states or 
maximally entangled states. Finally, we discuss the properties of MUBs consisting of 
product states only. We believe that understanding entanglement in MUBs can lead on 
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the practical side to novel applications and on the conceptual side to an understanding 
of why complete sets of MUBs can (not) exist for nonprime-power d. 

2. Conservation of entanglement 

Consider a bipartite system composed of subsystems A and B, i.e. its global dimension 
is d = dAds- Any (hypothetical) complete set of MUBs allows for efficient quantum 
tomography as it reveals complete information about an arbitrary quantum state of 
the system P, [15] . Hence we intuitively expect that the average entanglement over all 
the states constituting the complete set of MUBs shall be fixed with respect to some 
measure, independent of the choice of the bases. 

This intuition is made rigorous in this section. The relevant measure of 
entanglement is a function of the linear entropy of a reduced density operator. The idea 
of the proof is to use the property of a complete set of MUBs called a complex projective 
2-design [161 [17], which here means that the entanglement averaged over a complete set 
of MUBs is the same as the entanglement averaged over all pure quantum states. The 
latter is constant due to known results in statistical mechanics [18]. The message of this 
section, namely that the amount of entanglement is the same independent of a choice 
of the complete set of MUBs, may be well-known to scientists working with designs, 
but our proof is elementary and has immediate consequences for the distribution of 
entanglement among the states of MUBs. 

2.1. Complete sets of mutually unbiased bases and designs 

A complete set of MUBs is composed oi d+1 bases, each basis of d orthonormal vectors. 
We denote by the jth state of the mth basis, where for convenience we enumerate 
the states and the bases as j = 0, . . . , — 1 and m = 0, . . . , (i. To introduce the notion 
of a 2-design, one studies polynomials V{i) = P(xi, X2, i/i, 2/2^)5 which are biquadratic 
in variables Xi,X2 and separately in variables yl,y2, where Xi,yi are any coefficients of 
arbitrary state \i) with respect to a fixed (say, standard) basis and * denotes complex 
conjugation. Any complete set of MUBs is known to be a complex projective 2-design 
[T6l [TT] because the average of any V{jm) over states \jm) is the same as the average 
with the Haar measure over all pure states: 



2.2. The conservation law 

In order to utilize the design property of the complete set of MUBs in the studies of 
entanglement, we characterize the latter by the purity of a reduced density operator, 
say pA\jr^ = Trsdim) (jm|): 



(^(jm))MUBs = (^(^))h; 



(1) 



(2) 
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This quantity acquires its minimum of for maximally entangled states and its 
maximum of unity for unentangled product states. By 'maximally entangled states' 
we mean pure states with maximal possible entropy for the smaller of the subsystems. 
Note that due to the properties of the Schmidt decomposition it does not matter which 
subsystem is taken into account. Moreover, the assumptions behind Eq. ([T]) are fulfilled 
and, since per definition (P(j„))MUBs = ^(ipi) ^1=0 5]i=o Tr(p^l^.^^J, using the design 
property we write 

d d-l 

^ = E E ^<Pa\jJ = did + l)(Tr(pi|,))Haar. (3) 

m=0 j=0 

In the last step, we use the result by Lubkin ^18j, who studied how close the average 
reduced density operator is to a completely mixed state and found that 

(Tr(p^lJ)Haa. = (4) 

Therefore, the sum of entanglement over all the states of any complete set of MUBs is 
fixed and equal to 

£ = dAdB{dA + (Ib). (5) 

Note that the right-hand side is symmetric with respect to d^ and ds, which reflects 
the fact that we can as well study subsystem B. 

Eq. ([5]) has two immediate consequences. The first is that the distribution of 
entanglement among different states of a complete set of MUBs can be arbitrary as long 
as there is a proper amount of it. For example, Eq. ([5]) allows a complete set of MUBs 
to be formed by product and maximally entangled states as well as solely by partially 
entangled states. 

The second conclusion is that we cannot have a complete set of MUBs built entirely 
of product states or entirely of maximally entangled states. 

Lemma 1 Assume that dA ^ ds- In a complete set of MUBs which contains + 1 
product MUBs, all other bases contain only maximally entangled states. 

Proof. The sum of V{jm) over the states of product MUBs equals dAds^dA + !)• The 
only possibility to obtain the value of ([5]) is when for all the remaining d^dsids — 1) 
states, V{jm) acquires its minimal value of □ 

3. Complete sets of mutually unbiased bases in prime-squared dimension 

We showed that the complete set of MUBs may be chosen as consisting of product bases 
and bases containing only maximally entangled states. Here we present a construction 
of the complete sets with this property in dimension d = p"^ , where p is prime. 
The complete set will be generated from product MUBs with repeated application of 
a single entangling operation, in our case the control-phase operation. This makes 
our construction experiment-friendly. Explicit examples of MUBs generated by this 
method together with their factorization into product or maximally entangled bases are 
presented in the Appendices. 
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3.1. Complete sets of mutually unbiased bases in prime dimensions 

Before we present the new construction, let us briefly recall some of the known ones to 
which we will refer later on. If d = p is a prime number a complete set of p + 1 MUBs 
was first found by Ivanovic [2J. It is convenient to enumerate the bases as m = 0, . . . ,p 
with m = p corresponding to a standard basis, i.e. the basis in which the vectors of all 
other MUBs will be expressed. To simplify the notation and if no confusion arises, we 
will write the vectors of the standard basis without any index, i.e. |s) = \sp) enumerates 
the states of the standard basis. The other p MUBs have the Fourier-Gauss structure, 

1 ^"^ 

|j„^) = — for m = 0,...,p-l, and p > 2, (6) 

where ctp = exp(i27r/p) is the complex pth root of unity. The only exception to this 
formula is the case of p = 2 where one needs to refer to an imaginary unit i, the 
fourth rather than the square root of unity. For low dimensions, we present these bases 
explicitly in the Appendices. 

In odd-prime dimensions, a standard basis and a single MUB are sufficient to 
generate the complete set of MUBs with an application of a single unitary: 

W = diag[l,ap,al,...,a^;-'^% (7) 

which has the standard basis as the eigenbasis and permutes all other MUBs, i.e. 
W = \jm+i) with addition modulo p. 

Alternatively, one can construct complete sets of MUBs using Heisenberg-Weyl 
operators in prime dimensions, 

p— 1 p— 1 

X = J2\s + l){sl Z = J2c^»{sl (8) 

s=0 s=0 

with addition inside the kets modulo p. These operators span a unitary operator basis 
with respect to the trace scalar product as 

Tr[(X'^Z^)^X^Z'^] =p6a,cSb,ci. (9) 

According to the general result of Bandyopadhyay et al. [19] , if one can group elements of 
the unitary operator basis into disjoint subsets of d commuting operators (unity being 
the only common element of these sets), the common eigenbases of the commuting 
operators within each set are mutually unbiased. In the case of a system of a prime 
dimension, the groups of commuting operators can be chosen as powers of the operators 
Z, X, XZ, XZ^, . . ., XZP~^. Their eigenbases define a complete set of MUBs. It turns 
out that this set of MUBs is identical to the set of Eq. ([6]) up to the indexing of bases 
and states within bases. 

Lemma 2 Bases ([^ are the eigenbases of the operators X , XZ , . . ., XZ^~^ . 



Proof. Choosing the standard basis as the eigenbasis of Z, the eigenbasis of X is 
readily the Fourier basis, i.e. {|jo)}- Next note that for m = 1, . . . ,p — 1 we have 
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\jm) = ETJo ar™)^-''"^^ |.) with ^, = s + ... + {p-l) = l{p-s){p + s-l).The 
proof that these are exactly the eigenstates of XZ^™ is given in Ref. [19]. Since p is 
prime, 2m runs through all the powers of Z. □ 

These two methods of generating complete sets of MUBs in prime dimensions can 
be generalized to prime-power dimensions. However, these generalizations require a 
knowledge of elements of finite fields theory; see e.g. |3]. We now present our physically 
motivated construction of the complete set of MUBs in prime-squared dimensions d = p^. 

3.2. Two qubits 

We begin with a statement relating the number of MUBs to the possibility of swapping 
the states of subsystems. The statement itself holds for arbitrary dimension d = p'^, but 
it can be directly used to produce a complete set of MUBs of only two qubits. 

Lemma 3 Assume d = p^ and there exists unitary U that commutes with the swap 
operation S , and such that vectors {U \akhi)}, with k ^ I, form an MUB with respect to 
all product symmetric MUBs defined as {\amhm)}- Then {U \aibk)} is MUB with respect 
to all the bases mentioned above. 

Proof. The commutativity of U and S and the Hermiticity of S imply U = SU S. The 
assumed unbiasedness is expressed as |(am6m|f/|a'^6^)p = for all bases m = 0,...,p 
and all vectors \ambm) and \a'j^b'i). The computation of the overlap 

\{a,k\U^U\a[b[)\' = \{a,\a[){k\b',)\' = 1, (10) 

reveals that the basis from the thesis is unbiased to {U \akbi)}. The commutativity with 
the swap operation is used to prove its unbiasedness with respect to all product bases: 

\{ambm\U \a'ib'f.) ^ = \{ambm\SUS\a'ib'i,)\^ = \{am(3m\U\a'k(3'i)\^ = \, (11) 

where the last equality follows from the assumed unbiasedness and we put a = b, /3 = 
a, a' = b', (3' = a'. □ 

Note that the two bases {U \akbi)} and {U \aibk)} are simply related by the swap 
operation because U |a;&fc) = SUS laib^) = SU la^A) with a = b and f3 = a. 

In case of c? = 4, this lemma allows us to generate the complete set of MUBs 
starting with product MUBs. There are three MUBs in dimension 2 and therefore we 
begin with the following three product MUBs in dimension 4: {|ao&o)}, {^i^i)} and 
{10262)}- Consider now application of the control-phase (control-Z) operation 

V2 = ^{I ^ I + I ®a, + (T^^ I - a,^a,), (12) 

where I denotes a single qubit identity operator and ~ ^ q '^l ^ ^PP^y 

the control-phase operation on the two qubits prepared in states of the form |ao&i). 
The effect is best explained using Pauli operators. For a single qubit, we choose, in 
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accordance with Appendix A, the basis m = as the eigenbasis of = 

and basis m = 1 as the eigenbasis of = ^ ^ ^ ^ Therefore, the basis {|ao&i)} is 

the eigenbasis of the commuting operators (j^® I and I ^ ay, and their products. The 
control-phase operation maps these operators onto 

p2K®/)^2 = cr^®cT„ (13) 
V2{I ® (ry)V2 = a,^ ay. (14) 

The common eigenstates of these new operators are maximally entangled Bell states. 
Moreover, such Bell basis is mutually unbiased with respect to all our product MUBs. 
This can be verified directly or by using, e.g., the result of Bandyopadhyay et al. 
[T9] . We apply this theorem to tensor products of Pauli operators, eigenbases of 
which define our MUBs, i.e. the three product MUBs {|am&m)} are defined by sets 
of commuting operators {/ ® I, ax ® 1,1 ® (JxyCx ® {I ^ lyO'y (X> /,/ (8> (Jyjay ®ay} and 
{/ (g) /,(J2 ® /,/ ® cr^,cr2 cr^}, respectively, whereas the Bell basis V2 |ao^i) is defined by 
{/ ® I, ax ® az,az ® cry, ay ^ax}- Each set of four is clearly a set of commuting operators, 
and according to the mentioned theorem their eigenbases form MUBs. Since the V2 
operation is manifestly invariant under a swap of qubits, according to Lemma[3]we obtain 
the following complete set of MUBs: {|ao^o)}) {Wibi)}, {\(i2b2)}, {7^2 Wobi)}, {V2 |ai&o)}) 
which is explicitly presented in Appendix C. Note that for this dimension application 
of Lemma [3] has the same effect as the result of Ref. [21] stating for general dimension 
that if there is a set of d MUBs, then there also exists a set of + 1 of them. 



3.3. Two qupits 

Now we move to a system of a global dimension d = with p > 2. Two systems, each 
of prime dimension p, admit altogether + 1 MUBs. We shall show that they all can 
be generated via the multiple apphcation of a single entangling operation on product 
bases. For this purpose, we present a lemma which reduces the number of unbiasedness 
conditions one needs to check. 

Lemma 4 For p > 2 assume there exists unitary U such that \ {ambm\ U"' |ao&'„) P = 
for all < a,b,a' ,b' < p — 1, n = 1, ...,p — 1 and m = 0, ...,p and that [U, W I] = 
[U,I W] = 0, where W is defined in Eq. Then the bases {U" |a^6^+^)}, with 

fj,,!^ = 0,...,p — 1, together with the standard basis {\apbp)} form a complete set of 
MUBs. Addition of indices is modulo p. 

Proof. Consider an overlap between states of two bases of the proclaimed form 

M ^ I {a,b,U iUn^U-^' K'&W.') P = I U"'-" p. (15) 

Since U commutes with individual cycling unitary W ^ I and / ® W, it also commutes 
with their products. In particular, we have U"''" = {W^"' (g) ty'^'+^)f/'''-'"(W^-^' ® 
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y^-mu-uy We insert this expression into f llSp and since none of the bases there is the 
standard basis, the effect is to shift the indices of the local bases and get 

M = I {a.^.'b.^^'lU"'''' |a[,&;,_,> I' = ^, (16) 

where the last equality follows from our assumptions. Similarly, overlap with the 
standard basis equals | {apbp\ l^^^^^+y) P = ^ which follows from our assumptions 
after noting that the standard basis is not shifted by W , whereas the index of the other 
local bases we shift by — yU. □ 

Now we prove that the control-phase operation can be used to generate a complete 
set of MUBs in all prime-squared dimensions. The control-phase reads 

1 

V =-y a-'^^'Z'^ ® Z\ (17) 

^ a,fe=0 

Theorem 1 In every dimension d = with p > 2, there exists an integer 9 such that 
Vp satisfies requirements of Lemma \^ 

Proof. First note that since both Vp and W are diagonal in the standard basis, 
[pe^ W®I] = [V^,I®W] = is fuimiled for any 6. 

To prove that the bases {Vp"' \aQbn)} are unbiased to bases {\ambm)}, we refer 
once more to the results of Bandyopadhyay et al. [19j. They show that MUBs in 
prime dimensions may be chosen as eigenstates of sets of commuting operators 

X^Z^^" with /9 = 0, . . . ,p — 1 (see also Lemma [2]). The idea of the present proof is to 
show that operators defining bases {\aobn)} are transformed under the application of 
the control-phase into a new set of distinct operators which are all different from the 
operators defining bases {\ambm)}- Since commutation relations are preserved under 
unitary transformations, the results of [19] imply that the new operators define MUBs 
with respect to {\ambm)}- 

The control-phase acts symmetrically on both subsystems and we have up to a 
global phase: 

Since for different values of a and (3 the initial operators X° ® X^Z^^" were orthogonal 
with respect to the trace scalar product, the final operators are also orthogonal, i.e. 
we generated a set of trace-orthogonal operators which can be partitioned into proper 
commuting subsets. We now have to ensure that the generated set does not contain any 
operators determining product MUBs {\ambm)}- Since in the product MUBs the bases 
of A and B are the same, their defining feature is that operators determining the basis 
of A commute with the operators determining the basis of B. We check whether this 
commutation condition is satisfied by the operators on the right-hand side of (JTSll . The 
operators of A, i.e. X°Z^^", commute with the operators of B, i.e. x'^Z^'^""*""^'^, if and 
only if [19]: 

n{a'^9 + 2a(3 - fi'^O) = mod p. (19) 
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We are interested only in positive n and therefore ask whether the bracket is a multiple 
of a prime p. In other words, we are looking for the solution of the quadratic equation 
in the prime field J-p. It is well known that for p > 2 such equations have solutions if 
and only if there exists a field element 

A = 2(3Vl + e^. (20) 

Therefore, we need to choose such a value of 6 that 1 + 6"^ has no square root in the 
prime field J^p. Since for any element x in the field J-p we have = (p — x)^, we have 
no more than elements with square roots. That is, there exist an element x having 
a square root such that the next element, 1 + x, does not have a square root. Hence 
there always exists 9 such that \/l + 9"^ ^ J-'d- □ 

Generally speaking, there is no universal choice of 9 that is independent of p. We 
have neither found a function that for any given p returns 9 such that 1 + 9"^ does 
not have a square root in the field, and our construction generates the complete set of 
MUBs. A good guess of a useful value of 9 is often 1. Out of the first 1000 odd prime 
numbers, the construction with 9 = 1 fails in 494 cases, while out of the first 10000 odd 
primes it fails in 4988 cases. The lowest numbers for which this value does not produce 
the complete set of MUBs are 7, 17 and 23. 

3.4- Three qubits 

A similar construction using multiple application of only one entangling operation does 
not seem to exist for more than two subsystems of prime dimensionality. However, more 
operations can be used for the task. Here we show that three entangling gates can be 
used to produce a complete set of nine MUBs for three qubits. 

We start with the global standard basis {\abc)} and eight other bases that do not 
involve any local standard basis, i.e. \akbiCm) with k,l,m = 0, 1. We next apply to the 
basis {\akbiCm)} operation 

Qkim = ]^{I ® I ® I + Z'' ® ® Z"" + z^-^ ® z^-^ ® z^-™ -Z®Z®Z). 

(21) 

The resulting complete set of MUBs is given in Appendix F. 
3.5. Wocjan-Beth construction 

We would also like to mention the Wocjan-Beth construction [20], which is so far the 
only known construction that gives more MUBs in composite dimensions than there 
are for the smallest prime-power subsystem. The construction is designed for systems 
divisible into two identical subsystems. 

The method utilizes two kinds of vectors. The first kind is the so-called incident 
vectors, V. Exactly d of their d^ entries are equal to 1; the rest is 0. The task is to find 
families of d such vectors that satisfy the following requirements: within each family 
every pair of vectors is orthogonal, and two vectors from two different families have the 
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scalar product equal to one. For example, for d = 2 there are only three families of 
incident vectors, 
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} ■ (22) 



There is a one-to-one correspondence between families of incident vectors and mutually 
orthogonal Latin squares of order d. 

The other type of vectors is phase vectors, h, which have d complex entries, each 
of modulo 1. The two kinds are combined through operation "'|"". h 'I' V shall be 
understood as V with the first non-zero element multiplied by the first entry of h, the 
next by the second, etc. One needs d orthogonal vectors h and combines every phase 
vector with the incident vector using t- After normahzation we get as many MUBs as 
the number of incident vector families we found. 

We would like to mention that when we choose vectors h proportional to the rows of 
the Fourier matrix and the first two incident families in the most natural way (similarly 
to the example), two bases generated in this way possess a product structure whereas 
all others are maximally entangled. The present work suggests that it might be possible 
to extend this set with the 'missing' product bases, which would make the Wocjan-Beth 
construction even more powerful. 

4. Product mutually unbiased bases 

Our last topic is limitations on the number of MUBs and their entanglement, which 
follow from the fact that some bases are formed by product states. We call such bases 
product MUBs. First we present a straightforward bound on the maximal number 
of product MUBs; next we discuss classes of product MUBs to show that in every 
dimension one has two product MUBs such that there is no other product MUB with 
respect to them. There could still be entangled MUBs and we give an example in which 
this entanglement does not help us to build a complete set of MUBs. 

4.1. Maximal number 

We begin by showing that the only way to construct product MUBs in composite 
dimension dAds is to build them from MUBs in dimensions dA and ds separately. 

Lemma 5 Two product bases {\ab)} and {\a'b')} in dimension dAds are mutually 
unbiased if and only if \a) is mutually unbiased to \a') in dimension dA and \b) is 
mutually unbiased to \b') in dimension ds- 

Proof. If local bases are mutually unbiased, then clearly their product bases are 
also mutually unbiased. Conversely, assume the product bases are MUBs, i.e. 
|(a|a')p|(6|6')p = for all a,a',b,b' . Since the right-hand side is positive, neither of 
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the scalar products of the left-hand side is zero. In particular, this implies that keeping 
a, a', 6 fixed we have for all values of 6' that \{b\b')\'^ = l/dAdB\{a\a')\'^. Since the squared 
moduli are the quantum probabilities they sum up to = ^5 which implies 

that |(a|a')|2 = ^ and hence \ {b\b')\^ = ^- ^ 

The maximal number of product MUBs follows as a corollary. In a general 
dimension d = di . . .dn there are at most min^ Aij product MUBs, where Aij is the 
maximal number of MUBs in dimension dj. Note that the maximal number of product 
MUBs is also a corollary to Lemma [H 

4-2. Direct and indirect bases 

Not every set of product bases can be of the cardinality described below Lemma El The 
crucial distinction between the product bases is whether their states can be distinguished 
with (i) local measurements only or with (ii) additional classical communication [22] . 

The bases (i) are of form {\a) \b)} having all states \a) orthogonal in the first 
subspace and all states \b) orthogonal in the second subspace. We shall call them direct 
product bases because the matrix Pdirect having vectors \a) \b) as columns is a tensor 
product of matrices A and B having as columns vectors \a) and \b), respectively: 

^direct = A^B. (23) 

An example of a direct product basis is the tensor product of standard bases. 

Product bases (ii) can be written as \a) \b{a)), i.e. for every fixed vector \a) 
orthogonality of the product basis requires states of the second subsystem \b{a)) to be 
orthogonal, but importantly for different states \a) the orthogonal bases of the second 
subsystem may be different. The measurement in such product basis requires classical 
communication: after measuring the first subsystem the result needs to be fed-forward 
to a device measuring second subsystem in order to adapt its setting to a suitable basis. 
We shall call such bases indirect product bases. In matrix notation, the matrix of an 
indirect product basis Pmdirect cannot be written as a tensor product of matrices of local 
bases, but rather is of the form 

i^indirect = J]] |a) (a | O 5 (o) , (24) 

a 

where the columns of matrix B{a) are vectors \b{a)). An example of an indirect product 
basis of two qubits is 

|Oo)|Oo), |lo)|Oi), 

|Oo)|lo), |lo)|li). (25) 
4-3. Blocking product mutually unbiased bases 

A set of product MUBs is blocked if there exists no other mutually unbiased product 
basis with respect to this set. Indirect product bases lead to the minimal blocked set of 
product MUBs, and they have consequences for completeness of sets containing them. 
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Lemma 6 In every composite dimension, there is a blocked set of two product MUBs. 

Proof. The first basis is a standard basis: direct product basis. The second basis is 
an indirect product basis exhausting all possible MUBs for at least one subsystem. We 
order the subsystems such that di > d2 > ■ ■ ■ > dn- In dimension d2, there are at most 
(^2 MUBs with respect to the local standard basis. Since di > d2, for every orthogonal 
vector in dimension di one can have different orthogonal bases in dimension d2 which 
exhaust the whole set of local MUBs. According to Lemma |5] there is no other product 
MUB. □ 

A compact explicit example of an indirect product basis that together with the 
standard basis forms blocking product MUBs can be given in dimension being a power 
of a prime d = p^, which is regarded as the dimension of a Hilbert space of a set 
of r elementary p-level systems. For every subsystem there are exactly p MUBs with 
respect to the local standard basis and so is the number of distinguishable local states. 
Therefore, the basis {|(ji)o) I(j2)ji) |(j3)ji) • • • exhausts all allowed MUBs for 

all but the first subsystem. Here we denoted by \{jn)m) the state of the nth elementary 
subsystem in the mth MUB. 

The indirect product bases can block extendibility not only of a set of product 
MUBs but also of MUBs in general with no restriction to product bases. For example, 
in dimension 4 the set of three MUBs composed of the standard basis, the indirect 
product basis of Eq. ( 125|) and the Fourier basis cannot be extended by any other MUB 
[231 EU. 

5. Conclusions 

We were studying aspects of entanglement in states of MUBs in composite dimensions. 
Independently of the way a global system is split into subsystems, there is no complete 
set of MUBs that does not contain entanglement. In contrast, practically all MUBs 
are entangled as the dimension of at least one of the subsystems grows to infinity. The 
higher the dimension of the total system, the smaller the ratio of the number of product 
MUBs to the cardinality of the complete set of MUBs. This cardinality is proportional 
to the total dimension d, whereas the largest number of product MUBs is of the order 
of the smallest prime-power factor of d. Therefore, the ratio is the highest if ci = is a 
square of a prime, and even in this case the cardinality of the complete set is a square 
of the cardinality of the product MUBs and the ratio vanishes in the limit — )■ oo. 

We showed that entanglement of states of any complete set of MUBs is fixed. This 
has consequences for the distribution of entanglement among the states of a complete 
set and might be a useful hint for a search of the complete sets or one of the ingredients 
to (dis)prove their existence. This conservation law holds true independent of a division 
into subsystems and therefore perhaps an argument could be found that there is a finite 
set of divisions under which the entanglement cannot simultaneously match the proper 
value. Another route to follow is to begin with a set of states with a proper amount of 
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entanglement and apply local operations and classical communication in order to search 
for a complete set of MUBs. 

We also considered practical implementations of complete sets of MUBs and showed 
that for two subsystems, each with the same prime number of orthogonal states, the 
complete set can be generated via the multiple application of a single entangling 
operation on product states. The outcomes of this construction together with other 
examples of MUBs are explicitly presented in the Appendices for low dimensions (see 
also Ref. |25]). 

Turning to possible experiments, we note that there are various avenues for 
implementing quantum states in higher dimensions. For photons these include 
multiports and spatial-mode superpositions |26l [281 ISll |30l [31] or Hermite- 
Gauss and Laguerre- Gauss modes, most notably orbital angular momentum states 
[321 [331 [311 [351 [36]. 
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Appendix 

We present here explicit examples of complete sets of MUBs and for composite 
dimensions we emphasize division into entangled and product states. The notation 
used is explained on the example of a qubit {d = 2). 



Appendix A. d=2 

The symbol \jm) denotes the jth vector of the mth MUB. The standard basis is either 
denoted with subscript d or has no subscript at all: 

Note that when we write a basis as a matrix, we can freely permute columns, since it 
only changes the order of the vectors in the basis. 









Entanglement in mutually unbiased bases 



14 



Appendix B. d=3 
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where ct^ = exp27r/d is the complex dth root of unity. 



Appendix C. d=4 

The bases of this Appendix present exphcitly the result of construction described in 
section 3.2 of the main text. 
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|Oo)|Oi) + z|lo)|li) 
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|Oo)|Oi)-z|lo)|li) 
|Oo)|li)-^|lo)|Oi) 



where the kets refer to MUBs for the two-level system (Appendix A) and V2 is the 
control-phase operation between two qubits defined in Eq. f|T2|) of the main text. 
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Appendix E. d=6 

In this dimension it is not known if there exist more than three MUBs. A possible choice 
of three is to take the products 
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where the kets in two-dimensional vectors refer to qubit MUBs (Appendix A) and the 
kets in three-dimensional vectors refer to qutrit MUBs (Appendix B). 
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Appendix F. d=8 

The bases of this Appendix present exphcitly the result of construction described in 
section 3.4 of the main text. 
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The bases of this Appendix present exphcitly the result of construction described in 
section 3.3 of the main text for d = 3^ = 9. 



5n 



3 



1 

































1 

































1 

































1 



















|0) 













1 















■1 




ID 


> 














1 
















|2) 



















1 

































1 






























1 


) 










/ 


1 


1 


1 


1 


1 


1 


1 


1 




1 ^ 






1 


as 


ai 


1 


as 


ai 


1 


as 


ai 






1 


«i 


"3 


1 


al 


as 


1 




as 






1 


1 


1 


as 


as 


as 


ai 




«i 






1 




«i 


as 


al 


1 


ai 


1 




"3 






1 


al 


as 


as 


1 


ai 


ai 


as 


1 






1 


1 


1 


al 


al 


ai 


as 


as 


as 






1 


as 


al 


al 


1 


as 


as 


ai 


1 




V 


1 


al 


as 


al 


as 


1 


as 


1 








/ 


1 


1 


1 


1 


1 


1 




1 




1 


1 






03 


«i 


1 


as 


aA 


1 




as 




1 




0:3 


1 


al 


as 


1 


<4 




as 




1 


ai 




as 


as 


as 


al 


ai ai 




1 




1 


1 




aj 


1 


as 


1 


ag 


al 




as 


ai 


1 




al 


as 


1 


1 


ai as 




as 




1 


ai 




as 


as 


as 


1 


1 


1 






ai 


ai 




"i 


1 


as 


as 


a^ 


1 




1 


as 


«i 


V 


al 


as 


1 


as 


1 


ai 




1 


ai 


Q!3 / 




|Oo) ^ 




' |Oo) 


|lo) 


hi 


|lo) 


|2o) ^ 




, |2o) 




|0i) 

111) 

|2i) 



Entanglement in mutually unbiased bases 



/ 1 


1 


1 


1 


1 


1 


1 


1 


1 




1 




"3 


1 




a? 
"3 


1 


a 5 


"3 




1 






1 


a2 


Q/3 


1 


0:3 


0:3 


0:3 


1 


1 


1 


as 


as 


as 




«3 


1 


"3 


1 


as 


1 


as 


"3 


as 


1 


«i 


«i 


as 


1 


1 


«i 


as 


ai 


al 


al 


0:3 


as 


as 


1 


1 


1 


as 


al 


1 


1 


as 


al 


ai 


1 


as 


\ as 


1 




1 


al 


as 


ai 


as 


1 



1 

3 



1 
3 



/ 


1 


1 


1 


1 


1 


1 1 1 


1 \ 


as 


al 


1 


as 


1 


1 as ai 


1 


as 


1 


«i 


as 


"i 


ai as 1 


a 


2 
3 


1 


1 


1 


as 


1 


as ai ai 


a 


2 

3 


1 


as 


al 


as 


ai 


1 ai 1 


as 


al 


as 


1 


1 


«i 


as as 1 


a 


2 
3 


1 


1 


1 


al 


ai 


ai as as 


as 


al 


1 


as 


as 


ai 


1 1 as 


a 


2 
3 




1 


ai 


as 


al 


as 


1 as 1 


a 


11 






0o)|09 


) + as lo 


)|l9)+as|2o)|29) 


> 








0o)|09 


) + «i|lo 


)|l9) 


+ |2o)|29) 










Oo)|09) + |lo)|] 


-9) + 


«i|2o)|29) 










lo)|09) +as|2o)|l9) 


+ «3|Oo)|29) 






-. < 




lo)|09) + 


«3|2o)|l9) 


+ |Oo)|29) 




> = 






lo)|09) + |lo)|29) + 


«i|Oo)|29) 










2o)|09) + 


tt3|Oo)|l9) 


+ «3|lo)|29) 










2o)|09) + 


«i|Oo)|i9) 


+ |lo)|29) 










2o)|09) + |lo>|09) + 


«3|lo)|29> 






/ 


1 


1 


1 


1 


1 


111 


1 \ 


al 


1 


as 


al 


1 


as ai 1 


as 


ol 


as 


1 


al 


as 


1 ai as 


1 


1 


1 


1 


as 


as 


as ai ai 


a 


2 

3 


1 


as 


al 


as 


ai 


1 ai 1 


as 


as 


1 


al 


al 


as 


1 1 ai 


as 


1 


1 


1 


al 


ai 


ai as as 


as 


as 


al 


1 


1 


as 


ai ai 1 


as 




1 


al 


as 


ai 


as 


1 as 1 


a 


1/ 



Entanglement in mutually unbiased bases 



1 

71 



|Oo) 
|Oo) 
|Oo) 
|lo) 
|lo) 
|lo) 
|2o> 
|2o) 
|2o) 



09) + «i|lo)|29)+a^|2o)|l9) \ 



09)+a3|lo)|29) + 
09) + |lo)|29)+a3 



09)+«i|2o)|29)+a^|0o)|l9) 



09)+a3|2o)|29) + 
09) + |2o)|29)+a3 



09)+ai|0o)|29)+Q;i|lo>|l9 



09)+a3|0o)|29) + 
09) + |0o)|29)+a3 



2o) 
2o) 



Oo) 
Oo) 



lo) 
lo) 



I9) 
I9) 



I9) 
I9) 



I9) 
I9) 



/ 


1 


1 


1 


1 


1 


1 


1 


1 


1 \ 




1 


as 


al 


1 


"3 


al 


1 


as 


al 




1 


al 


as 


1 


al 


as 


1 


al 


as 




as 


as 


as 


al 


al 


al 


1 


1 


1 




al 


1 


as 


1 


as 


al 


as 


al 


1 




1 


al 


as 


as 


1 


al 


al 


as 


1 




as 


as 


as 


1 


1 


1 


al 


al 


«i 




1 


as 


al 


al 


1 


«3 


as 


al 


1 


V 


al 


as 


1 


as 


1 


al 


1 




^3 / 



1 

71 



|0i)|09) + |ll)|29) + 

|0i)|09)+a^|li)|29) 
|0i)|09) +a3|li)|29) 
|li)|09) + |2i)|29) + 
|li)|09)+al|2i)|29) 
|li)|09)+a3|2i)|29) 
|2i)|09) + |0i)|29) + 
|2i)|09)+ai|0i)|29) 
|2i)|09)+a3|0i)|29) 



|2l)|l9) 

+ a3|2i) 
+ «i|2i) 
|0i)|l9) 
+ a3|0i) 
+ al\0,) 

\h)\h) 

+ Ois\ll) 

+ al\h) 



I9) 
I9) 

I9) 
I9) 

I9) 
I9) 



Vs 



/ 


1 


1 


1 


1 


1 


1 


1 


1 


1 \ 




al 


1 


as 


al 


1 


as 


al 


1 


0:3 




al 


as 


1 


al 


as 


1 


al 


as 


1 




as 


as 


as 


al 


al 


al 


1 


1 


1 




al 


1 


as 


1 


as 


al 


as 


"i 


1 




as 


1 


al 


al 


as 


1 


1 


al 


as 




as 


as 


as 


1 


1 


1 


"i 


al 


al 




as 


al 


1 


1 


"3 


al 


"3 


1 


as 


\ 


al 


as 


1 


as 


1 


al 


1 


«i 


"3 / 



Entanglement in mutually unbiased bases 



1 

71 



|0i) 
|0i) 
|0i) 

111) 
111) 
111) 

|2i) 
|2i) 
|2i) 



09) + ai|li)|l9)+«^|2i)|29) 



09) + |ll)|l9)+«3 
09)+a3|ll)|l9) + 



09)+ai|2i)|l9)+a^|0i)|29) 
09) + |2i)|l9)+a3|0i) 
09)+a3|2i)|l9) + |0i) 
O9) +Q;i|0i)|l9) +« 
09) + |0l)|l9)+a3 
09)+Q;3|0i)|l9) + 



2i) 
2i) 



4\ 

ii; 
ii; 



29) 
29) 



29) 
29) 

ll)|29) 

29) 
29) 



1 
71 



IO2) 
IO2) 
IO2) 
II2) 
II2) 
II2) 
122) 
122) 
122) 



09) + 112)119) + 122)129) 
l9) + ai|22) 
I9) +a3|22) 
) + 102)129) 
h) + al\02) 
I9) +a3|02) 
> + 112)129) 

l9)+«i|l2) 
I9) +Q;3|l2) 



O9) +a3|l2) 

09)+«i|l2) 
09) + |22)|lc 
09)+«3|22) 

09)+ ail 22) 
09) + 1 02)119 

09)+a3|02) 
09)+«i|02) 



29) 
29) 

29) 
29) 

29) 
29) 



/ 1 


1 


1 


1 


1 


1 


1 


1 


1 


\ 


1 


as 


a\ 


1 


as 


"3 


1 


as 


ai 




1 


al 


^3 


1 


al 


as 


1 


ai 


as 




"i 


al 




1 


1 


1 


as 


as 


as 




as 


al 


1 


al 


1 


as 


1 


as 


ai 




1 


al 




as 


1 


ai 


ai 


as 


1 




«i 


al 




as 


as 


as 


1 


1 


1 




1 


"3 


«3 


al 


1 


as 


as 


ai 


1 




\ "3 


1 


a\ 


1 


al 


as 


ai 


as 


1 


/ 



1 1 


1 


1 


1 


1 


1 


1 


1 


1 


\ 


as 


ai 


1 


as 


ai 


1 


as 


ai 


1 




as 


1 


ai 


as 


1 


ai 


as 


1 


ai 




ai 


ai 


ai 


1 


1 


1 


as 


as 


as 




as 


ai 


1 


ai 


1 


as 


1 


as 


ai 




ai 


as 


1 


1 


ai 


as 


as 


1 


«i 




ai 


ai 


ai 


as 


as 


as 


1 


1 


1 




ai 


1 


as 


as 


ai 


1 


1 


as 


«i 




\ as 


1 


ai 


1 


ai 


as 


ai 


as 


1 


) 



Entanglement in mutually unbiased bases 



24 



1 

7! 



^2/ 


IOq^ 


1 


"0 1 -"-2/ 1 -^9/ 


+ 2o) lo) 
-p Lj; j 1 ^2/ 1 -'-9/ 






1 


\U)\2q) + 

1 -"-2/ 1^9/ ~ 


"3 1^2/ 1 -"-9/ 


2/ 


IOq^ 
1 9/ 


1 


"3 1 -^2/1^9/ 


+ |2o)I1q) 
~ \ -"-9/ 


I2) 


IO9) 

1 


+ 


"3122)129) 


+ as|02)|lq) 


I2) 


IO9) 


+ 


122)129) + 


«I|02)|l9) > 


I2) 


IO9) 


+ 


«i|22)|29) 


+102)119) 


22) 


IO9) 


+ 


«3|02)|29) 


+ a3|l2)|l9) 


22) 


IO9) 


+ 


102)129)+ 


«^|l2)|l9) 


22) 


IO9) 


+ 


«3|02)|29) 


+112)119) , 



P3 



02) 1 


[ |0i) 


12) 


5 111) 


22) J 


[ |2i) 



where the kets refer to MUBs for a qutrit (Appendix B), and Vj, is the control-phase 
gate for two qutrits as given in Eq. ( |T7I) of the main text. 
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